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ABSTRACT 



We study the reaction kinetics of end-functionalized polymer chains dispersed in an unre- 
active polymer melt. Starting from an infinite hierarchy of coupled equations for many-chain 
correlation functions, a closed equation is derived for the 2nd order rate constant k after pos- 
tulating simple physical bounds. Our results generalize previous 2-chain treatments (valid 
in dilute reactants limit) by Doi de Gennes and Friedman and O'Shaughnessy 0, 
to arbitrary initial reactive group density no and local chemical reactivity Q. Simple mean 
field (MF) kinetics apply at short times, k Q. For high Q, a transition occurs to diffusion- 
controlled (DC) kinetics with k x^/t (where Xt is rms monomer displacement in time t) 
leading to a density decay ~ no — n^xl. If no exceeds the chain overlap threshold, this 
behavior is followed by a regime where Ut ~ ^/x^ during which k has the same power law 
dependence in time, k ^ xf/t, but possibly different numerical coefficient. For unentangled 
melts this gives n^ ~ t^^/^ while for entangled cases one or more of the successive regimes 
Ut ~ t~^^'^, t~^/® and t~^^^ may be realized depending on the magnitudes of Q and no. Ki- 
netics at times longer than the longest polymer relaxation time r are always MF. If a DC 
regime has developed before r then the long time rate constant is /c ~ R^/t where R is the 
coil radius. We propose measuring the above kinetics in a model experiment where radical 
end groups are generated by photolysis. 



PACS numbers: 82.35. +t, 05.40.+j, 05.70.Ln 
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1. Introduction 



The study of polymer-polymer reaction kinetics is a fundamental problem in polymer 
science. Such reactions occur in many technologically important processes such as vulcan- 
ization of rubbers and free radical polymerization [Q. Compared to analogous reactions 
between small molecules, polymer reaction kinetics are novel in that they reflect static and 
dynamic properties of the polymer chains which are host to the reacting groups. Theoretical 
studies predict that rate constants depend both on the degree of polymerization and time 
. Model experiments and numerical simulations |]|,|] testing these laws thus offer 
a way to probe fundamentals of polymer dynamics . In this paper we focus on reactions be- 
tween end-functionalized chains in a polymer melt. So far, no systematic experiment exists 
measuring rate constants in a melt as a function of molecular weight and time. It is our aim 
to motivate such model experiments by developing a complete theoretical picture for such 
reactions for direct comparison with experiment. 

The situation we analyze is illustrated in fig. ^ We imagine that at t = a certain 
fraction of the chains randomly and uniformly distributed in a monodisperse polymer melt 
carry chemically reactive end-groups. Experimentally this can be realized for example by 
attaching photocleavable groups to the ends of a certain fraction of the melt chains |T0| . 
A laser pulse then cleaves these groups into radical pairs (fig. 0), one of which is a small 
molecule and the other of which is attached to the end of the chain (a macroradical) |TT|] . It 
was shown in ref. |T^ that after a transient the small radicals disappear, leaving behind a 



fraction of order unity of the macroradicals. The reaction kinetics of these kinetically isolated 
macroradicals will then follow the kinetics developed in the following. 

Anticipating second order kinetics, our aim is to determine the second order time-dependent 
rate constant kt defined by 

. _dnt ^ 
= ^ = -kt rit , (1) 

where rit is the number density of reactive groups at time t. Now for the usual case of 
reactions between small molecules, kt is independent of time. This is so because small 
molecules obey simple Fickian diffusion, Xt ~ t^^'^, where Xt is the rms displacement in time 
t. This is a very dilute exploration of space (the exploration volume x^ increases faster 
than t) and diffusion is fast enough to smooth out any reaction-induced correlations p|,|12|. 



Thus reactant distribution remains uniform and random at all times, the number of pairs 
per unit volume which are in contact at time t is therefore a^nf, and hence nt ~ Qa^ri^ 
where the local reactivity Q measures probability of reaction per unit time when in contact. 
This implies a time-independent rate constant k ^ Qa^. For "infinitely" reactive groups for 
which Q is effectively equal to the rate at which diffusion brings apart a pair of molecules 
in contact, i.e. Q ~ where ta is the reactant relaxation time, this leads to the well 
known Smoluchowski rate constant k ^ /ta ~ -Dsmaii^ where -Dgmaii ~ /ta is the 



small molecule self-diffusivity. 

Interpolymeric reaction kinetics are more complicated. For times shorter than the longest 
polymer relaxation time r, corresponding to a diffusion distance of order the coil radius i?, 
reactive groups attached to polymer chains explore space in a dense way: The dynamical 
exponent defined by Xt ~ t^^^ ^ is 4 or 8 depending on the degree of entanglement and time 
n,r3|; thus the exploration volume x^ is increasing more slowly than t and the number of 



times reactive groups visit points within their exploration volumes is an increasing function of 
time. This was clearly pointed out by de Gennes P[ who differentiated between "compact'' 
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exploration when z is greater than the dimensionahty of space d and "non-compact" for 
z < d. 

Now in the hmit in which reactants are very dilute, n^^B? <^ 1, the probability that by 
time r a reactive group's exploration volume contains a second reactive group is very small. 
Hence for t < t reactions are due to the few isolated pairs which happened to be initially 
within diffusive range, i. e. the members of which were initially within Xt of one another. Due 
to the compact exploration, for times greater than t the number of collisions between such 
a pair grows without bound. Therefore for "infinitely" reactive groups {Qta ~ 1) reaction is 
almost certain by t. Thus a depletion hole of size Xt starts to grow in the 2-body correlation 
function. The number of reacted pairs per unit volume is x^n^, which when differentiated 
with respect to time implies a short time time- dependent rate constant kt ~ a;^/t. Diffusion 
dynamics change at r, since center of gravity Fickian diffusion of the whole chain takes over, 
Xt ~ t^/^. The hole in the 2-body correlation function stops growing and the situation is 
as in the small molecule infinite reactivity case, but now with an effective capture radius R 
reflecting the short time compact dynamics. Thus 

f x^/t {t < t) 

kt ~ \ {Qta = 1, n^R" < 1) . (2) 

[ R'/r it > r) 

The t > T expression for k was flrst obtained by Doi [jl6|,|l| for unentangled melts and was 



extended by de Gennes to t < r and entangled cases. Friedman and O'Shaughnessy 
derived k from first principles using renormalization group calculations. 

The generalization of eq. (Q) to the concentrated regime, no-R^ ^ 1, is far less obvious 
and is one of the main objectives of this paper. In this case the timescale by which a reactive 
group will almost certainly contain another group within its exploration volume occurs during 
the compact regime, t < t. During this regime one can no longer assign reactions to 
pairs since a given group has the opportunity to react with more than one partner. A 
proper study of the reaction kinetics requires a full many-body treatment. Such a many- 
body formalism involves an infinite hierarchy of coupled dynamical equations for correlation 
functions of different orders [|T7| , p!8| . Other than renormalization group analyses |ll9| , p0| . 



this complication is typically resolved by decoupling approximations expressing higher order 



correlations in terms of lower order correlation functions pT|. In our approach we are able 
to derive closed equations after assuming much less restrictive bounds on the magnitude of 
correlation functions. We show that the implicit assumption of ref. ^ namely that eq. @ 
is also valid in the concentrated regime, is correct. 

Our study is also general in terms of the local reactivity Q. This is very important 
experimentally, since typical chemical reactivities are extremely small ||2^j23|| , Qta ~ 10~^, 
with the exception of radical species which are very close to being infinitely reactive, Qta ~ 1 



In the following section we present the many-body formalism describing polymer-polymer 
reaction kinetics and introduce our assumptions on the magnitude of correlation functions. 
This results in a closed equation for k. We solve this equation in section 3 in terms of Xt 
and identify sequences of kinetic regimes depending on the values of Q and no- We apply 
these general results to the case of unentangled and entangled melts in sections 4 and 5, 
respectively. Results are presented in the form of a "phase-diagram" in the Q-uq plane, 
different regions of which correspond to different kinetic sequences. We conclude with a 
discussion of our results in section 6. 
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2. Solution for Rate Constant kt 



Consider tlie situation illustrated in fig. |I| in which reactive ends are initially randomly 
distributed with density tiq. Reactions commence at t = 0. Let us define ptiji, to be the 
2-body correlation function of reactive chain ends located at ri, r2. Translational invariance 
implies pt is a function of ri-r2 only. Now the reaction rate per unit volume at point r 
is proportional to the number density of pairs in contact at r, a^p(r,r), multiplied by the 
local reactivity Q. Here a is the monomer size. Using the fact that translational invariance 
implies the average density nt is independent of r one has 

rit = -Xpt{0,0) , X = Qa^ . (3) 

Thus the calculation of the reaction rate requires evaluation of the 2-body correlation func- 
tion. A rigorous dynamical equation for pt in the case of reactions between small molecules 
has been derived by Doi in ref. In the polymer reaction problem considered here, in order 
to derive a closed relationship for pt in terms of the degrees of freedom specifying the location 
of the reactive ends only, one must first average out the degrees of freedom specifying the 
locations of the other monomers; this is non-trivial and requires renormalization group (RG) 
methods. However, RG studies of 2-body bulk polymer reaction kinetics l2^J^ indicate that 



the basic physics is completely captured by the approximate closing of the system in terms 
of the coordinates of reactive groups: correct scaling behaviors are obtained, only the pref- 
actors being unreliable. Making the approximation of closure in terms of reactive degrees 
of freedom, then the pt dynamics involves the 3-body reactive groups correlation function 
(ri,r2,r3) MM- 



Pt(ri,r2) =nl-X J dr[dr'^ dt' Gt^t>{ri,r2,r[,r'^) ptir[,r2)6{r[ - r'^) 

- X j dY\ dr'2 c/r'g J dt' Gt-t'{ri, r2, r'^, r'^) p'^\r[, r'^, r'^) {6{r[ - r'g) + S^r'^ - r'^)} , 

(4) 

where Gi(ri, r2, r'j^, is the equilibrium chain end propagator giving the net weighting for 
two ends to arrive at ri,r2 given starting points ^[,^'2 in the absence of reactions. This is 
a well known object. The sink terms on the right hand side of eq. (4) (illustrated in fig. 
^ describe the three ways in which reactions diminish pt from its initial value tiq. The first 
two-body sink term subtracts off pairs which failed to reach ri,r2 because their members 
reacted with one another at r'^ at time t'. The remaining two sink terms subtract off pairs 
which would be at ri, r2 but only one member of which reacted at time t' at location rg. Such 
a reaction involves a third chain, weighted by the appropriate 3-chain correlation function. 
These are many- chain terms; were they absent, one would have a closed 2-chain system. In 
refs. D and ^ these terms were omitted. As explained in the introduction one expects such a 
2-body approach to be vahd in the dilute limit. In eq. (4) we used a 6 function as a reactive 
sink. This is a coarse-grained description of the reaction process over a scale of order the 
monomer size a. 

Notice that eq. (4) is not closed in terms of pt since it involves the unknown p\ . It is 
in fact impossible to write a closed exact equation for pt since correlation functions of all 



orders are coupled in an infinite hierarchy of dynamical equations [|18| , [r7| . This complication 
which arises in all many-body reacting systems is typically resolved by approximating 3- 
body correlations in terms of lower order correlation functions In this study we are able 
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to write a closed equation for pt after assuming less restrictive bounds on the magnitude of 
correlation functions. To explore the constraints imposed by these bounds we first transform 
eq. (4) to an integral expression for the function gt(r,0) = pt{r,0)/nf. Notice that eq. (|I]) 
then implies kt = Agt(0, 0). It is demonstrated in appendix A that after making this change 
of variables, one has 

qtir, 0) = 1 - I dr[ dt'Gt-t'ir, 0, r'^, v[) kt> 

+ 2 j rfr;rfr'2|^*rft'Gt_i,(r,0,r; + r^,r'2)A:i,{pHr;|0)-pf (r;|0,0)} , 

(5) 

where pt(r|0) and pp^(r|0, 0) are the conditional reactive chain end densities at r, given one 
and two reactive ends, respectively, at the origin: 

p,WO).^^, pf)wo.o).^i!i^ . (6) 

Eq. (5) can be reexpressed in terms of the propagator G'^'^^(ri, = / dr^Gtiv'-^^Tc'^ + 
ri, r'g, + r2), namely the probability density two chains ends are separated by ri at time t 
given initial separation V2- Thus eq. (5) becomes 

gi(r,0) = 1 - fdt'GlZ'ir,0) k, + ^^(r) , (7) 

where 

^,(r) = f^dt' j dv'G]Z'{v,v')k,pt'{r) , ^,(r) = 2 [pt{v\Q) - pf\v\Q,Q)] . 



We want to solve eq. (|^) for kt or equivalently for ^^(0,0). This requires information on 
the properties of v?t(0) which in turn involves 2-body and 3-body conditional densities. Now 
we make the assumption that the more reactive groups placed at the origin, the lower the 
conditional density. Chemical reactivity can only induce anticorrelations. Thus: 

pP''(r|0,0) < Pt(r|0) <nt (assumption) . (9) 

Eq. @ immediately implies the following constraint on p (see eq. (8)) 

< ptij) < 2nt . (10) 

Eq. (P) also implies that < 1. Thus the magnitude of the second term on the rhs of eq. 
(1^) must exceed that of the 3rd term, dt'G^^]!^,{r,Q)kt' > V't(r). This is the full inequality 
which must be obeyed by the function pt' (r') involving the unknown p^^^ . Integrating over r 
yields the following less demanding inequality: 

t dt'kf > t dt' j dr' kf pt'ir') . (11) 
Jo Jo J 
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On the strength of the above constraints, eqs. (p!0| ) and (|Tl|), we will now argue that the 
solution of eq. (|^) for gj(0,0), but with the term (pt{0) deleted, gives the correct power law 
solution for kt to within a constant prefactor. Expressing qt{0, 0) in terms of kt, setting r = 
and deleting v?t(0), eq. (0) becomes 

kt = \-\ f dt'St-t'kt> , (12) 

JO 

where we have introduced the return probability St = Gt^^{0, 0). 

In appendix B, a self-consistent argument is presented to justify the deletion of (pt{0). A 
summary of this argument is as follows. Now if one accepts eq. (0), then eqs. (|T^) and 
(|l|) lead to a sequence of power law regimes in time, both for kt and rit (these are explicitly 
obtained in the next sections). Using these solutions in the constraints (eqs. (|ToD and (|Tl|)) 
and in the expression for v't(O) (eq. (8)), the function /i™^^(r') which maximizes (pt{0) (for 
a given time t) subject to the constraints is determined. This in turn implies an upper 
bound on v^t(O), namely (y9™^^(0), according to eq. (8). It is then shown in appendix B 
that (f'^^^i^O) -C 1. (Of course, in reality there is one unique value for (pt{0), since /it'(r') 
is a determined function, eq. (8). What has been achieved here is to bound this function 
given our incomplete knowledge.) Thus ift{0) may be deleted in eq. (|^ without error; even 
prefactors are expected to be correct. There is, however, one exception to this simple state of 
affairs. During those time regimes where rit ~ 1/xt (see next sections) we find ip^^^{0) = A, 
where A is a constant of order unity. For these regimes, we are forced to make one further 
assumption, that kt remains a power law during this regime. In appendix B, it is shown that 
this implies that the possibility of ift{0) being of order unity may modify the coefficient of 
kt only, not the power itself. 



Thus in the following we simply solve eq. (|T^)- This equation has been derived in ref. 



starting from a 2-body formalism valid in the dilute reactive species limit as discussed. Here 
we argue that it is valid for all uq. Expressions equivalent to eq. (O) were the starting point 



in ref. ^ the authors of which pointed out the existence of a rit ^ regime in polymer 



reaction kinetics. 



Now Laplace transforming t ^ E one immediately solves eq. (|T^) for k{E): 



*(^' = WTxsm ' '''' 



which is the main result of this section. 



3. Kinetic Regimes and Timescales 

In the previous section we obtained a solution for k in Laplace space. This involves the 
return probability St- It is easy to show that St ~ ^/xt where Xt is rms displacement. 
In this section by substitution in eq. (|13D we obtain solutions for k and n in terms of Xt- 
Explicit forms of Xt for unentangled and entangled melts are considered in sections 4 and 5. 

For times shorter than the longest polymer relaxation time r (corresponding to E ^ t), 
Xt ~ t^^^ with z = 4 and 8 depending on time and the degree of entanglement P^ , P^ . Thus 
St ~ t~^^^ which implies S{E) ~ ii^^/z-i -^jj^ ^ ^_ Yqt times t ^ t, center of gravity 
Fickian diffusion applies, Xt ~ t^^^- Thus St ~ t^^^^ and it follows that for 1/E ^ t, S{E) 
approaches its = limit, dtS(t) ~ t/R^ (see refs. |l],|I). Thus, as a function of l/E, 
S{E) increases and then saturates at r. 
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To determine kt we must determine which of the two terms, 1 or XS{E) in eq. (0) is 
dominant. There are two cases. If Q (or A) is sufficiently small, Q < Q (see definition 
below), then XS{E) -C 1 for all E values, implying 

kt^X, {Q<Q) . (14) 

For Q > Q there exists a 1/E value (shorter than r), corresponding to a timescale after 
which XS{E) > 1. In this case during the regime < E'^ < r, k{E) ^ l/ES{E) ~ E-^/\ 
Thus Laplace inverting eq. (|T3|) we have: 

{t < t*) 

(t^<t<r) , (Q>Q) • (15) 

{t » r) 

The quantities ^2 and Q are defined by 

^ = A (t; < r) , Q = ^ . (16) 




The value of ^2 is determined by demanding continuity of kt- In Laplace space, corresponds 
to the 1/E value value at which XS{E) becomes of order unity. Eq. (|TBp has a solution for 
^2 only if Q is large enough so that ^ < r. Thus the condition t2 = T defines Q in eq. (|TB]). 

Now let us see what eqs. (0) and (p!5D imply for the decay of rit, considering the cases 
Q > Q and Q < Q in turn. 

1. Case Q > Q. The solution of eq. (0) is rit = no/{l + rio /q kfdt'). Thus using eq. (|T5|) 
one has 

l/{l + t/t*J {t<^t*) 

1/(1 + noX?) (t^«t«r) iQ>Q) . (17) 



no 



1/(1 + VCi?) 



We define the timescales t*^, t*^ ^ and by 



1 



- Ano ' ^"^'^ " noi?3 ' ~ no ' ^^^^ 

Here ti is the time to diffuse a distance of the order of the typical initial separation between 
reactants. Notice that at a certain timescale the time-dependent term in the denominator of 
rit in eq. ([THp becomes of order unity. This time is one of t*^,ti, and t^i?- Different kinetic 
behavior occurs depending on the ordering of these timescales amongst themselves and with 
respect to ^2 and r. 

Let us consider first the kinetics in the "concentrated" case, no-R^ > 1, i.e. U > t. Then 
the monotonically increasing time-dependent term in the denominator of rit in eq. (p^Tj) 
becomes of order unity before r. There are two possibilities depending on whether or not Q 
is large enough such that is smaller than ti (notice that ti is independent of Q). Now if 
indeed < U one has 

^ no - An^t ^ no - n^x^ ^ 1 ^ [t; < U < r, S=°°^) . (19) 

We use the symbol 8*^°'^'^ for "strong concentrated" since eq. (0) applies for high Q and no 
values (see phase-diagrams of figs. H). Notice that we do not need to specify Q > Q m. eq. 
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(|T9|) since only for such Q values is < t- If on the other hand ^ > ti, one similarly has 

nt^no- Xnlt 1 ^ 1 -1. ^J-j- (Q > Q, ti < f,, W) (20) 

where W stands for "weak" since these systems are less reactive when compared to the 
strong concentrated regime (see figs. ^ and |^). Notice that we do not need to specify ti > t 
in eq. (^), since this is implied by Q > Q (or equivalently < r) and ti < 

Now in the the "dilute" case, uqR^ < 1, the time-dependent term in the denominator of 
fit in eq. ([17D becomes of order unity for t > r: 



t* . t t* l ^ 

nt ^ no -Xnlt ^0-^00;^^ nQ~nlR^- {Q > Q, k > t, S<^'^) . 

r K'^t/T 

(21) 

Here S'^'' stands for "strong dilute." 

2. Case Q < Q. Finally let us examine the Q < Q case. Substituting eq. (0) in (|ID 
one immediately has: 

nt ^ no- Xnlt ^ ^ {Q<Q, WW) . (22) 

At 

We use the symbol WW (very weak) to refer to systems belonging to this class which 
correspond to the lowest Q values in which mean field kinetics are always applicable. 



We have thus identified 4 distinct behaviors, eqs. (p!9|)-(|22|), for the decay of nt, corre- 
sponding to 4 regions in the Q-no plane. In the following two sections we explicitly plot 
these regions for both unentangled and entangled melts after giving expressions for xt as 
dictated by the Rouse and reptation model, respectively. 

4. Application to Unentangled Melts 

The dynamics of polymers in an unentangled melt, i. e. with degree of polymerization 
shorter than the entanglement threshold A'^e; are well known to obey Rouse dynamics |jl5|,|l^ 




-=N\ - = N'/' . (23) 

ta a 

Thus Rouse dynamics is characterized by compact {z = 4) short time behavior, followed by 
noncompact Fickian diffusion {z = 2). 

By substitution of eq. (^3|) in eqs. ( |16|) and ( |l8l) we obtain explicit expressions for ^2 and 

tf. 

t; = ta/{Qtaf iQ>Q) , ti = tJia'no)^/' (noR' > 1) (24) 

(the timescale ti is only relevant to the kinetics when no-R^ > 1). Timescales t"^, t*^^ are 
already given explicitly in eq. (p^S]). Using these expressions, in fig. ^jthe lines ^2 = t, ti = r, 
and ti = t^ have been drawn in the Q-Uq plane. Certain sections of these lines are omitted, 
in those regions where they are irrelevant. These lines define four distinct regions. It is 
straightforward to verify that each one of these corresponds to one of the kinetic regimes 
defined in eqs. (|l^)-(0) as indicated in the figure. 
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Notice that the condition ti = ^ also imphes ti = ^ = t^, as one can easily verify using 
eqs. (0) and (|TBp. In the Q-uq plane, the = = tm ^^^^ separating the S™'^'^ from the W 
region thus defines a density-dependent reactivity Q*: 

Q*t, = Ka3)i/3 ^ (25) 

indicated in fig. ^ 

Perhaps the most interesting feature of the unentangled case is the rit ~ t^^/^ decay of rit 
during the rif ^ ^/xf regime in regions 8'^°'^'^, W. 



5. Application to Entangled Melts 



In the case of very long chains, N > N^, polymer motion is affected by entanglements. In 



the reptation model p!5|JT^ entanglements are assumed to inhibit lateral chain motion on the 
scale of the "tube" of diameter re = N^^'^a, corresponding to a portion of chain comprising 
Ne units. For times shorter than the diffusion time te = N^ta to distance r^, monomers do 
not feel the tube and obey Rouse-like t^^^ dynamics as in unentangled melts. For t > te, 
the chain diffuses curvilinearly up and down the tube in 1-dimensional t^^^ Rouse motion. 
During these "breathing modes," monomer rms displacement in space increases as smce 
the tube is itself a random walk. The chain relaxes its configuration relative to the tube by 
time tb = N'^ia, corresponding to monomer diffusion distance r-^ = {N/Nef^. For Ion ger 
times, coherent diffusion along the tube gives rise to t^^^ rms monomer displacement. This 
regime persists until the longest polymer relaxation or "reptation" time, r = {R/r\y)H\y = 
{N^ /Ne)ta, by which time the chain has completely diffused out of its initial tube into a 
new and uncorrelated one (here R = N^/'^a). The process then repeats itself indefinitely, 
corresponding to long time Fickian center of gravity motion, Xt = R^t/rY^"^. In summary. 



a [t/taf^ {t<te = NX) 

re {t/tcY^^ {te<t<h = NX) 

rb (t/tb)'/' {t^<t<T = (R/r^)^ = NX/Ne) 

R {t/T)^l^ it > t) 



(26) 



Thus, there are 3 compact regimes with a sequence of dynamical exponents z 
followed by noncompact z = 2 exploration. 
Substituting eq. (EBI) in eq. (P^) one obtains 



4,8,4, 



taiQta)-' 

te{Qtea^/rl)-^"> 
%{Qhayrl)-^ 



{Q > r!/{tea')) 

(rl/Xa^) <Q< r^e/itea^)) 

{Q<Q< rl/Xa^)) 



(27) 



Notice that, unlike the unentangled case, the appropriate formula for ^2 is now dependent 
on the dynamical regime during which it happens to occur. This is determined by Q. Similar 
remarks apply to t;, expressions for which depend on hq. Eq. (^61) in eq. (llSD leads to 



tainoa')'/' 
te(norf)^/2 



(no > ^) 



(r^ ^ >no>r^^) 



(28) 



> no> R~ 
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Once again the condition ti = ^ = t*^ defines a density-dependent reactivity Q*\ 

[ {noa^y/yta (no > r-^) 

Q* = \ (nor3)8/3/(t,noa3) {r-^>no>r^') (29) 
[ {noriy/y{t^noa^) (r^ ' > uq > R'^) 

Similarly to the unentangled case, the phase diagram of fig. ^ is constructed. Now the 
three Q > Q regions develop fine structure. In each of these subregions, crossovers in rit 
behavior occur during different short time compact regimes. 

How does one use this phase diagram to determine the reaction kinetics for a given system? 
The system's Q and no values define a point in the diagram. Depending on which of the 4 
regions, S™"'^,S'^'', W, WW, this point happens to belong to, the reaction kinetics are then 
given by the appropriate member of the 4 equations (p!9|)-(p2D. But in these expressions 
for reaction kinetics timescales ^ind ti may appear. The appropriate formulae for these 
timescales are determined by the fine structure location of the point within the region (eqs. 
( p^ and (p^)). The other timescales, tm'^m.R' if relevant, are given by eq. (|T8]). 

For example, let us consider point a in fig. ^belonging to region 8'^°'^'^. According to eq. 
(|19D the relevant timescales (in addition to r) are ^ and ti. Fig. ^ then indicates < and 
te < ti < ti,. Thus the Hq — Xn^t regime of eq. (|1^) occurs before t^, followed by Uq — ri^x^ 
which occurs during the latter part of the first Xt ~ t^^^ regime and the earlier part of of the 
Xt ~ t^^^ regime. The crossover to rit ~ ^/x^ occurs during the t^^^ regime, lasting up to r. 
During these times rit ~ t^^^^ followed by rit ~ t~^^'^. 



6. Conclusions 



This study has addressed the kinetics of reactions between end-functionalized polymer 
chains diffusing in an unreactive polymer melt matrix, as a function of initial reactant density 
Uq and chemical reactivity Q. Results were presented in the form of a "phase diagram" in 
the Q-uq plane. 

At short times, simple mean field (MF) kinetics apply: the rate constant k ^ Qa"^ is 
independent of time. These kinetics give rise to a density decay nt ~ no — rio t/t*^. In regions 
W and WW of the Q-rio plane (figs. ^ and ^ these MF kinetics persist for times greater 
than t*^, leading to nt ~ l/Qt. Physically, as defined in eq. (0) is the timescale after 
which a reactive group is very likely to have reacted with the "mean reactive field" supplied 
by all other reactants in the system. 

A transition to diffusion-controlled (DC) kinetics occurs for sufficiently reactive groups, 
Q > Q. This is a result of compact exploration of space at times shorter than the longest 
chain relaxation time r. Thus a timescale exists after which any pair which was initially 
within diffusive range (i. e. the exploration volumes of its members are overlapping by t) is 
bound to react. If occurs at times shorter than the time ti needed for a reactive end to 
diffuse a distance of order the typical separation between reactive groups (regions S"^"^^, S'^'^ 
of figs. ^ and §), then the depletion in the initial density no is proportional to the number 
of pairs with initial separation less than xf. rit ^ no — no^^. Equivalently, kt ~ Xt/t p[. 
For unentangled melts this leads to kt ~ t~^^^, while for entangled cases successive regimes 
kt ~ t~^^^, and t~^^^ may exist depending on Q, uq. 

In regions S*^°"'^ and W kinetics are diffusion-controlled for times longer than ti. During 
these times the exploration volumes of reactive groups overlap. We find a density decay 
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Ut ~ l/^^t 10 • Roughly, this means that only one reactant exists within a region of size 
comparable to the exploration volume x^; had there been several, they would have reacted. 
For unentangled melts this implies Ut ~ while an extra Ut ~ t~^/^ regime arises in 

entangled systems. These DC kinetics are analogous to the anomalous long time decay 
Ut ~ l/a^t in the case of small molecule reactions (A + A — > 0) obeying Fickian dynamics, 
Xt ~ t^^^, in spatial dimensions below a critical dimension dc = 2. For arbitrary dynamics, 
Xt ~ t^^^, the critical dimension is = z; thus c? = 3 is below dc for t < r since z is then 
either 4 or 8. Notice that the rate constant during these times, kf ^ xf/t, is the same as 
for short time {t < ti) DC kinetics (as described above). However there is no obvious reason 
to expect them to be identical since each one has different physical origins. Indeed, our 
calculations suggest that the numerical coefficients of these two rate constants are different. 

All diffusion-controlled kinetics are truncated at r which marks a crossover to Fickian 
center of gravity polymer diffusion. Thus long time kinetics are always MF. In cases where 
DC kinetics preceded, the short time kinetics are reflected in a "renormalized" rate constant 
k ^ R^/t |]T|,0] leading to rit ^ l/{Rh/T) at long times. 

Experimentally, the most interesting region of the Q-hq plane is 8*^°'^'^. This region in- 
cludes the rather peculiar nt ~ ^/xf decay. Probing reaction kinetics in this region would 
require extremely reactive species such as radicals (or electronically excited groups as in 
photophysical systems 0) since typical chemical reactivities, Qta ~ 10~^, are far below the 
Q = Q threshold, even for very long chains (see figs. ^ and fig. One can imagine exper- 
imentally probing the 8^°"'^ region by photocleaving functional groups attached to polymer 
chain ends with a laser pulse, generating radical pairs with initial density no as in fig |^. 
Each pair consists of a polymer radical (P) and a monomeric radical (M). An experiment 
of this type was studied theoretically in refs. [1^ and There are 2 underlying principles: 
(1) For large N, since polymers are much less mobile than monomers, the polymer-polymer 
reaction constant kpp is much less than the monomer-monomer constant /cmm- (2) Polymer- 
monomer reactions are dominated by the more mobile monomer [^, i.e. /cmp ~ ^mm to 
within a prefactor of order unity. In the following we take kMP = kMM for simplicity. It was 
shown in refs. that after a time T ^ {kMAino)^^ virtually all M have reacted, leaving 

behind a finite fraction of the P: the polymer radicals are kinetically isolated and the t > T 
kinetics involve polymer-polymer reactions only (fig. [l|). Specifically, after time T a fraction 
1/e of the P remain whilst only a fraction e/e of the M remain, where e = kpp /kMM ^ 1- 

Returning to the present case, if we assume all radical groups have similar reactivities 
Q ^ we identify T ^ (Ario)"^ = ^m- (This is the MF timescale we have met previously, 
as expected: small molecules in d = 3 obey MF kinetics.) Let us suppose no is chosen so 
the polymer radicals are strongly overlapping, uqR^ ^ 1 as in fig. ^ The conclusion is that 
after our reactive polymer system belongs to the 8'^°'^'^ region [^] since Q exceeds the Q* 
threshold. But in region S^°'^'^ by definition is much less than the time ti for a reactive 
polymer end group to diffuse the distance between reactive ends. Thus (consistent with the 
above arguments) macroradicals cannot have substantially reacted with one another by time 
t^. Monitoring the decay of macroradical density for t > t^, the full sequence of kinetics is 
thus experimentally accessible, i. e. rit ^ uq — n^xl for t*^ < t < ti followed by n^ ^ for 
ti < t < T. Experiments of this type thus hold promise for the measurement of fundamental 
scaling laws of polymer-polymer reaction kinetics. 

This work was supported by the National Science Foundation under grant no. DMR- 
9403566. 
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Appendix A. Derivation of eq. (5) 



Eq. (4) can be written equivalently as |T8 



^^^^dt^"^^ -^Pt(ri,r2) = -Api(ri,r2)(5(ri-r2)-A j dv^pf\vi,Y2,V3) {6{ri - rg) + 6{y2 - rg)} 

(Al) 

where the terms on the right hand side represent reaction sinks and V is the diffusion 
operator. The propagator G in eq. (4) is the inverse of (d/dt — V). Notice that —Vp is equal 
to the divergence of the reactive chain end current j at ri, Vp = — (Vn + Vr2)j(i"i, i"2)- 
(For example for small molecules one would have j = — D(Vri + ^r2)P) where D is the 
diffusion coefficient.) Since j must be linear in p, one has Vpt = T>qtn'f = n^Vqt. Then, 
substituting pt{ri,r2) = qt{ri,r2)n'f in eq. ( |A1| ) one obtains 



%(ri,r2) 
dt 



Vqt{ri,V2) = -Agt(ri,r2)(5(ri - ra) + 2pt(ri - r2, 0)ktnt 



2p?\r, 



r2,0,0)/n2 
(A2) 



where we have performed the integration in eq. ([Al| ) and used the identities Pti^i, r2) = 
Pt(ri — r2, 0) and pf\ri, r2, ri) = p[^\ri, r2, r2) = pP^(ri — r2, 0, 0) which result from trans- 
lational invariance. Inverting eq. (A2) one obtains eq. (5) after using eqs. (|l|), and the 
definitions of the conditional densities in eq. (0). 



Appendix B. Bounds on ^pti^) 



This appendix shows that the deletion of the positive term ft{0), which results in eq. (|^ 
resulting in eq. (IT) gives the correct scaling solution for kf. We show this by first determining 
an upper bound, namely (y9™^'^(0), by finding the function /i™'^^(r') in the definition of if ( eq. 
(8)), which subject to the constraints satisfied by the real p, eqs. (|TD|) and (|TT]), maximizes 
(pt{0)- We consider each of the 4 regions in figs. |^, |^ in turn. 



1. Region S'^™'^ 

A. Times t <^ti. During these times rit ~ uq to within higher order corrections (see eq. 
(p!9D). In this case we just need to consider the constraint imposed on p by eq. ( p!OD only 
(including constraint ( [Tl| ) would lead to an even more strict bound which is not necessary 
here). Thus eq. (|T0|) implies p^^'^iv') = 2nt. Substituting ^"^.x ^g-j performing 
the r' integration leads to 

after use of eq. (p!5|). In the S™'^'^ regime <^ and t ^ U. Thus eq. ( PT| ) implies 
^max^g^ -C 1. Heucc duriug this regime (y? is unimportant since it is much smaller than terms 
1 and g((0,0) in eq. (|^). Eq. (|T^) is thus expected to give the correct k (even the prefactor 
will not be modified by deletion of ip). 

B. Times U <^ t <^ t. In this case we consider both constraints ([T0|) and (pA)). Now 
clearly from its definition in eq. (8), V'f(O) is maximized when the integral of k{t')pti{r') 
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has its maximum value; this corresponds to replacing the inequality in (^TJ) by an equal- 
ity. Moreover, /ij'(r') should be distributed around the points at which Gl'^J^^,{0,r') in the 
definition of v?f(0) in eq. (8) is maximum. But G^'^^ has the well known scaling form 




GnO,r)^K/"* (B2) 



Thus G'j!_''j/(0, r') is roughly constant (in space) for r' < Xt-t' and vanishes for larger r' . 
For a given r', it achieves maximum values at the maximum available t' values outside the 
vanishing region. Furthermore, we note that the maximum amplitude of /i™'^^ is of order 
n{t'). Hence the function maximizing V5t(0) is approximately 

/^r^(r') ^ I ^^'^ ® ~ "'^ ^ , jl k{t')dt' = J' dt' J dr'A:(t')/xr^(r') . 

(B3) 

where G is the step function (0(x) = for x < O,0(x) = 1 for x > 0). That is, /i™^^(r') 
is localized at (t', r') close to (t, 0) where it has its maximum amplitude. The time tc 
defines the lower limit of the support of ^"^^^ and is determined by demanding that the 
total integral of k^^^^ yields the required value. Assuming that tc is of order t (which is 
verified self-consistently), we determine tc by substituting kt and rit from eqs. (|15D and 
([191) in eq. (|B^ ) and performing the r' integration. It is simple to show that this implies 
jt^jti^ — uf'l^ judu ~ 1, after using Xi ~ t^l^ . Hence tc = Pt, where /5 is a constant of order 
unity (0 < /? < 1). 

Substituting /x™'^^ into the expression for v^t(O) in eq. (8) and integrating over r' one 
obtains ^^'■''{0) ^ Jl dt' k{t')n{t') ^ \n[n{tc) / n{t)\ after use of eq. (|). Since tc = (3t thus 
^max^g^ _ where A is a constant of order unity. 

We now argue that, assuming kt is still a power law, and given the boundndess of V5t(0), 
this implies kt ~ x^/t during this regime. Substituting kt ~ in the integral term in eq. 
(|^) and setting r = 0, one easily derives the time-dependence of this term is ^""(3-^)/^. Since 
(^t(O) positive and bounded by A, and since the term gf(0,0) tends to at long times, the 
only way to satisfy eq. (0) for r = is by setting a = {?> — z)/ z. This implies kt ~ x^/t. For 
this value of a the integral term in eq. (|^) tends to a constant. Hence if </5t(0) is of order 
unity during this regime, it must be a constant in order for eq. (^ to be satisfied. The 
actual numerical prefactor of kt will then be a function of the value of v^f(O). 

C. Times t ^ t. During these times rit ~ l/{B?t/T) and kt ~ B? /t (see eqs. (|19|) and 
(|15|)). We use the same reasoning as in part l.B. up to eq. (|B3|) . As in I.B., it is easy to 
show that tc/t = 1 — 7, where 7 ^ [i?^t/(a;fr)]^/^ and that v^™^(0) ^ 7. Since 7 <^ 1 for 
t ^ r, it follows that ^tiS^) ^ 1 and that V5t(0) may be deleted from eq. (|^. 



2. Region W 

A. Times t <^t*^. Similar arguments as in part l.A. apply. 

B. Times t*^ <^ t <^ t"^. The procedure is the exactly the same as in part l.C, but now 
using kt ^ X and rit ~ 1/At (see eqs. (p!5|) and (^0]) ). In this case one has tc/t = 1 — 7 with 
7 ^ {Xt/x^Y^^ and ipf'^^^i^O) ^ 7. Since 7 <^ 1 for t <^ t2 it follows that v^t(O) ^ 1 during 
this regime. 

C. Times t > tg- Identical to l.B. for t^ '^t <^t and to part l.C. for t:^r. 
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3. Region S^^'^ 



A. Times t j^. For times t <^t, same as part l.A. During this regime <y9™'^^(0) in 
eq. ( [BID is much smaller than unity since ti > r. For r ^ t ^ t^^, similarly to l.A., one 
finds ip^'^^{{]) = t/t*^^ which is much smaller than unity. 

B. Times t^t*^^. Same as part l.C. 

4. Region WW 

A. Times t <^ t*^. Same as part l.A. 

B. Times t*^. Same as part 2.B., but now notice that 7 -C 1 for all t {t^ is not defined 
in this region). 
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Figures 

FIG. 1. End-functionalized polymer chains (coil radius R) dispersed within a polymer melt 
of inert but otherwise identical chains. The initial density of reactive groups is no. 

FIG. 2. Schematic of the situation immediately after photolysis of functional end-groups by a 
laser pulse generating radical pairs (consisting of a monomeric radical and a macroradical) . 
The concentrated regime (overlapping reactive coils) is illustrated. After a transient most 
of the more mobile monomeric radicals react leaving behind a fraction of order unity of 
the polymer radicals. The situation of fig. |I] is then recovered. 

FIG. 3. The depletion in the number density of reactive groups at ri, r2 originates from two 
terms in eq. (4). The two-body term subtracts off those pairs which would have been at 
ri,r2 at time t, but failed to do so because both members reacted at an earlier time t' at 
point The many-body term subtracts off pairs which would have been at ri,r2 had 
there been no reactions, but failed to arrive because one member of the pair reacted at 
an earlier time. 

FIG. 4. (a) Unentangled melts, reaction kinetics "phase diagram" in Q-Uq plane. Axes 
are logarithmic and units chosen such that = a = 1. Maximum possible density is 
nott^ = (all chains functionalized) . Different sequences of kinetic behavior arise in 
the 4 regions S^°^^, S'^'', W and WW. (b) As (a), but reactivities and densities expressed 
in terms of degree of polymerization N. 

FIG. 5. As fig. ^, but for entangled melts. Regions S™"^"^, S'^'^ and W now develop fine 
structure. In a given sub-region each relevant timescale occurs within a given reptation 
diffusion regime thus defining a unique sequence of kinetic regimes. 
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Fig. 3 
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Fig. 4(b) 
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Fig. 5(b) 



